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ABSTRACT 

> A  modular  method  for  calculating  hydrodynamic  interactions  between 
particles  in  low-Reynolds-number  flow  has  been  constructed  by  using  multipole 
expansion  solutions  for  the  reflection  field.  The  approach  is  made  possible 
by  the  use  of  Faxen  laws  in  relating  the  multipole  moment  to  the  incident 
field.  The  method  is  illustrated  and  checked  by  recalculating  known 
expressions  for  the  resistance  and  mobility  tensors  for  two  spheres.  The 
method  can  be  readily  generalized  to  handle  three-particle  (or  n-particle) 
interactions  as  shown  in  a  following  paper.  New  forms  of  the  Faxen  laws  for 
prolate  spheroids  are  given  and  will  form  the  basis  for  other  papers  on 
spheroid-spheroid  and  spheroid-wall  hydrodynamic  interactions.  The  important 
result  is  that  "first-reflection*  solutions  can  be  readily  calculated  even  in 
cases  where  the  ambient  velocity  field  is  obtained  by  a  numerical  procedure. 
These  results,  as  asymptotic  (far-field)  solutions,  furnish  a  check  for  more 
robust  codes.  In  addition, -they  are  important  on  their  own  since  the^  provide 
crucial  information  for  the  renormalization  theory  used  in  suspension 
rheology.  ^ — _ _ _ 
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SIGNIFICANCE  AND  EXPLANATION 


The  calculation  of  hydrodynamic  interactions  between  particles  is  needed 
for  the  understanding  and  control  of  many  natural  and  manufacturing  processes, 
for  instance,  those  involving  sedimentation,  colloidal  stability  or  suspension 
rheology.  In  these  applications  the  external  forces,  torques  and  dipole 
moments  on  the  particles  are  known  a  priori  and  the  problem  is  to  calculate 
the  resulting  translational  and  rotational  motions.  In  practice,  since  the 
governing  differential  equation  requires  knowledge  of  these  motions  for  the 
boundary  condition,  one  has  to  solve  first  for  the  forces,  torques  and  dipole 
moment  in  a  collection  of  translational  and  rotational  problems  and  then 
invert  to  obtain  the  desired  motions. 

The  purpose  of  this  paper  is  to  show  that  these  problems  can  be  solved 
directly.  Explicit  calculations  and  comparisons  with  other  techniques  are 
shown  for  interactions  between  spheres.  The  first  step  towards  the 
corresponding  calculation  for  spheroids  (which  arise  in  the  modeling  of 
suspension  of  nonspherical  particles,  e.g.  clay/water  systems)  are  also  given. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


HYDRODYNAMIC  INTERACTIONS  BETWEEN  PARTICLES 
IN  LOW -REYNOLDS -NUMBER  FLOW: 

A  MODULAR  APPROACH 


Sangtae  Kim 

1.  INTRODUCTION 

Hydrodynamic  interactions  (HI)  appear  in  the  analysis  of  many  different  problems  such 
as  kinetic  theory  of  polymers/  mechanics  of  suspensions,  and  the  convection  of  particles  in 
porous  media.  In  this  series  of  papers,  a  modular  approach  to  hydrodynamic  interactions 
between  particles  of  arbitrary  shapes  in  the  creeping  flow  regime  (vanishing  Reynolds 
number)  is  introduced  and  applied  to  problems  in  each  of  these  areas.  The  approach  will  be 
illustrated  by  extending  the  literature  on  HI  in  two  directions:  interactions  involving 
axisymmetric  particles  (prolate  spheroids)  and  interactions  involving  Store  than  two 
particles  (three  spheres  in  a  uniform  stream). 

The  governing  equations  for  the  fluid  (viscosity  u)  velocity,  v  and  pressure,  p 
are  taken  as 

(1.1)  -*p  ♦  U?2v  -  0  , 

with  boundary  conditions  at  the  surface  of  the  particle 

v  •  y  ♦  «  *  4 

where  U  and  w  are  the  particle  translational  and  rotational  velocities.  The  geometry 
is  either  that  of  a  number  of  particles  in  a  fluid  of  infinite  extent  as  shown  in  Figure  1, 
a  single  particle  in  a  domain  with  boundaries  as  shown  in  Figure  2,  or  a  combination  of  the 
two. 

The  approach  is  modular  in  the  sense  that  the  HI  solution  is  obtained  by  a  two  step 
procedure  —  a  calculation  of  properties  intrinsic  to  each  interacting  particle,  followed 
by  a  method  for  integrating  these  modules.  The  power  of  the  method  is  due  to  the  fact  that 
once  the  first  step  has  been  accomplished,  the  second  step  follows  as  a  'back  of  the 
envelope*  calculation.  As  illustrated  in  later  sections,  this  feature  leads  to 
considerable  savings  in  human  computation,  particularly  for  the  'stobillty  problems"  for  the 
calculation  of  the  motion  of  particles  under  imposed  forces  and  torques. 


Sponsored  by  the  United  States  Army  under  Contract  No.  r>AAG29-80K>004t. 


V 


The  method  is  a  generalization  of  the  HI  calculation  for  spheres  using  the  method  of 
reflections  as  given  in  Happel  and  Brenner  ( 1 °65 ) ,  and  thus  finds  similar  utility  and 
limitations.  The  limitations  are  imposed  by  the  complexity  of  the  algebraic  manipulations 
for  the  higher  order  contributions.  Therefore,  accurate  solutions  to  problems  involving 
strong  interactions  as  defined  by  Ganatos  et  al  (1980)  are  precluded.  However,  the  method 
gives  accurate  results  for  weak  interactions  and  thus  one  immediate  application  is  its 
utility  in  the  testing  of  more  robust  numerical  codes  in  this  limiting  case.  In  addition, 
there  are  fundamental  problems  In  suspension  mechanics  where  the  leading  order  terms  in  HI 
as  determined  from  the  method  of  reflections  plays  an  important  role  (for  examples  on 
suspensions  of  spheres,  see  Batchelor  1972,  Batchelor  and  Green  1972b,  and  Hinch  1977). 

The  analogous  calculations  for  particles  of  arbitrary  shape  can  be  accomplished  using  the 
techniques  presented  here. 

The  work  is  presented  in  four  parts,  starting  with  this  introductory  paper.  In  part 
II,  the  method  is  used  to  calculate  the  interaction  between  two  prolate  spheroids,  which  in 
turn  leads  to  a  generalization  of  the  Sotne-Prager-Yamakawa  potential  used  in  the  kinetic 
theory  of  polymers.  In  part  III,  calculations  for  the  hydrodynamic  interaction  between 
three  spheres  illustrate  the  utility  of  the  method  for  multi-particle  systems.  The  inter¬ 
action  between  a  convected  particle  and  the  solid  matrix  of  a  porous  medium  is  the  subject 
of  part  IV. 

The  following  section  is  an  exposition  on  the  method  of  reflections,  as  used  in  this 
series  of  papers.  The  multipole  expansion  solution  for  the  velocity  field  around  the 
particle  is  introduced.  A  procedure  for  calculating  these  multipole  moments,  for  particles 
of  arbitrary  shape,  is  shown  in  Section  3.  It  is  a  variation  of  the  procedure  for  deriving 
the  "Faxen  laws"  discussed  in  earlier  works  (Faxen  1922,  1927,  Brenner  1964b,  Batchelor  and 
Green  1972a,  Rallison  1978).  The  present  procedure  exploits  certain  reciprocity  relation¬ 
ships  between  the  Faxen  laws  and  velocity  representations  in  certain  associated  flow 
fields.  The  method  is  illustrated  by  deriving  new  and  possibly  more  useful  forms  of  the 
Faxen  laws  for  prolate  spheroids. 
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In  Section  4  and  5,  the  "module  integration"  idea  is  applied  to  various  problems 
involving  HI  between  spheres.  Most  of  the  results  presented  in  these  sections  are  already 


in  the  literature.  However,  the  calculations,  in  addition  to  providing  a  test  of  the 
method,  explain  the  surprising  occurrences  of  zero  coefficients  in  the  expansion  for  the 
mobility  tensor.  Furthermore,  the  mobility  results  are  obtained  directly,  without  the 
usual  steps  of  taking  linear  combinations  of  subsidiary  problems. 

2.  THE  METHOD  OF  REFLECTIONS 

The  basic  solution  strategy  is  similar  to  the  method  of  reflections  as  given  in  Rappel 
and  Brenner  (1965).  Their  *  and  **  notation  is  replaced  with  a  system  of  subscripts 
1  and  2  because  the  latter  are  more  readily  generalized  for  multi-sphere  interaction 
problems.  The  original  aspect  of  this  work  is  the  development  and  use  of  Faxen  laws  which 
greatly  facilitate  the  method. 

For  two  widely  separated  particles,  with  centers  of  reference  at  jg^  and  £2'  the 
zero-th  order  solution  for  the  velocity  field  is  simply  the  sum  of  the  disturbance 
solutions  for  each  particle  in  isolation,  i.e.  without  hydrodynamic  interaction.  In  the 
present  notation,  these  single-particle  solutions  are  denoted  v1  and  v2 .  Since  the  no¬ 
slip  boundary  conditions  on  each  particle  are  violated  by  the  velocity  field  emanating  from 
the  other  particle,  we  correct  the  solution  by  adding  on  two  new  fields  that  cancel  these 
discrepancies.  However,  any  field  that  helps  with  the  boundary  condition  at  one  particle 
will  upset  matters  at  the  other  particle,  hence  we  get  a  sequence  of  velocity  fields 
comprising  an  iterative  approximation. 

At  each  reflection,  the  fields  that  are  created  to  help  satisfy  the  boundary 
conditions  are  known  as  reflection  fields.  The  field  that  violated  the  boundary  condition 
(hence  forcing  the  creation  of  the  reflection  field)  is  called  the  incident  field.  A 
schematic  representation  of  the  reflection  process  is  given  in  Figure  3.  The  following 
convention  will  he  used  to  label  the  subscripts.  For  a  reflection  at  particle  B , 

(B  ”  t,  2)  the  reflected  field  will  be  labeled  by  adding  the  subscript  B  to  the 
subscripts  of  the  incident  field.  Note  that  the  isolated,  single-perticle  solutions  may  be 
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F1GUKE  3. 

Schematic  representation  of  the  reflection  process  for  two  spheres. 


considered  as  reflections  from  the  ambient  velocity  field  v  .  This  will  be  called  the 
zero-th  reflection  in  accord  with  the  literature.  The  fields  reflected  from  the  single- 
particle  solutions  are  the  first  reflections,  and  those  reflected  from  the  n-th  reflection 
fields  are  referred  to  as  the  (n+1)-th  reflection  fields. 

Given  an  incident  field,  we  need  a  method  for  calculating  the  reflected  field.  For 
spheres,  one  can  use  Lamb's  general  solution,  with  Hobson's  (1955)  addition  theorems  for 
transforming  the  spherical  harmonics  from  a  coordinate  centered  at  one  sphere  to  one 
centered  at  the  other,  as  shown  by  Happel  and  Brenner  (1965).  Instead,  we  present  an 
alternative  approach  which  can  be  applied  to  more  general  particle  shapes.  This  approach 
is  based  on  the  Faxen  laws,  the  integral  representation  for  the  Stokes  solutions,  and  the 
multipole  expansion. 

An  incident  field  v1 ( x)  which  disturbs  the  boundary  condition  at  a  sphere  is 
countered  with  a  reflection  field  having  the  following  integral  representation  (see  Howells 
1974): 

(2.1)  -  ^  ^(O-n)  •  I(x-x')  dS(x')  . 

Here  a  is  the  stress  evaluated  at  the  surface,  jc'  is  a  vector  from  the  center  of 
reference  to  a  point  on  the  surface,  and  I  is  the  Green's  dyadic  or  Oseen  tensor  given  by 


(2.2) 


W** "  M  4ij +  Vj  • 

The  multipole  expansion  solution  is  obtained  by  expanding  the  Oseen  tensor  about  the  center 
of  reference.  Once  the  multipole  moments,  £  (o»n)x’ ...x*  dS(x'),  are  related  to  v1, 
the  solution  can  be  carried  out  by  calculating  each  reflection  to  the  desired  accuracy. 

Such  relations,  known  as  Faxen  laws  are  available  for  spheres  (Faxen  1922,  1927, 
Batchelor  and  Green  1972)  and  for  ellipsoids  (Brenner  1964b,  Rallison  1978).  Therefore, 
the  reflection  procedure  just  described,  can  be  applied  readily  to  particles  of  these 
shapes.  For  nonspherical  particles,  however,  resolution  problems  can  arise  in  the  module- 
integration  step.  To  improve  the  resolution  for  nonspherical  particles,  alternate  forms 
for  the  Faxen  laws  will  be  derived  in  the  following  section. 
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3.  CONSTRUCTION  OF  THE  FAXEN  LAWS 


From  the  previous  section,  it  follows  that  the  reflection  fields  can  be  constructed  if 
the  multipole  moments  can  be  expressed  in  terms  of  the  incident  field.  This  section 
presents  the  Faxen  laws  for  the  monopole  and  dipole  moments,  i.e.  the  force,  torque  and 
stresslet,  on  a  particle  of  arhitrary  shape  movinq  with  a  translational  velocity  and  an 
angular  velocity  jo.  Since  the  idea  for  this  method  came  from  the  striking  similarity 
between  the  Faxen  laws  and  Stokes  solutions  for  a  sphere,  the  method  will  be  introduced  by 
applying  it  to  spheres. 

The  Faxen  laws  for  the  force  (Faxen  1922),  torque  (Faxen  1927)  and  stresslet 

00 

(Batchelor  and  Green  1972a)  on  a  sphere  in  an  ambient  field  v  are  compared  below  with 

«D 

the  velocity  fields  for  a  stationary  sphere  in  a  uniform  stream  U  t  vorticity  field  « 
and  rate  of  strain  field  5  (Hinch  1977): 

Force: 

a^  2  “ 

(3.1a)  F  =  6*tia(1  +  7  )v  (x^  -  6*uaU  . 

Solution  for  uniform  stream: 


(3.1b) 


v(x) 


-  6rija0  •  ( 1  +  V2  ) 

O 


Torque : 

(3.2a) 

Solution  for  vorticity  field 
(3,2b) 

Stresslet : 

(3.3a) 


T  =  4xuaJ(7xv  (x^)  -  2u) 


t  3  t 

v(x)  =  ft  #x  +  4wua  ft  *V 


8*ii 


i-r*-3"  +tov2)  i  (V^  +  ^Mx-x, 


Solution  for  rate  of  strain  field: 
(3.3b)  v(x) 


20  3  a2  2  i  (~”~1 1 

Ev*  +  _  spa  (K-V)-(l  ♦  w  V  ) 


This  similarity,  as  noted  by  Hinch  (1977)  is  rooted  in  the  Lorentz  reciprocal  theorem  for 
Stokes  flow.  In  fact,  the  reciprocal  theorem  can  be  used  to  generate  the  Faxen  laws  as 
shown  by  Brenner  (1964b)  and  Fallison  (1978).  However,  this  formal  procedure  can  be 
streamlined  by  using  the  fundamental  solutions  for  the  velocity  as  shown  below. 


The  single-sphere  aolution  (3.1b)  satisfies  tha  boundary  condition  at  r  •  a: 


(3.4) 


6»ga(  1  +  ’2) 


§5u 


If  we  go  back  to  the  integral  repreaentation  for  a  aphere  in  an  ambient  velocity  field 
v"  (x) 


(3.5)  v(jt)  -  v“(x)  -  —■  £g  (O'n)'l(x-JS' )  dS(x') 

and  then  operate  on  both  aidea  with  the  LHS  of  (3.4),  the  result  la 

2  2  m 

(3.6)  6wga(  1  ♦  7  IjjtlS,)  -  6»ga(1  ♦  V  )v  (2l)  -  o*^  dS(x' )  . 

In  the  aphere,  Z  ~  51  *  )  *°  that 

2 

6»ua(1  +  ~~  ^2)v(x,)  “  6*gaU  . 

The  reeult  ia  (3.1a),  the  Faxen  law  for  the  force  on  the  aphere.  Looking  back,  we  aee  that 
the  procedure  is  equivalent  to  Brenner's  (1964b)  approach  since  the  integral  representation 
is  a  special  case  of  the  reciprocal  theorem. 

The  tensor  identities  which  follow  from  the  boundary  conditiona  for  (3.2b)  and  (3.3b) 
similarly  give  the  Faxen  lawa  for  the  torque  and  streaalet  when  applied  to  the  integral 
repr eaentat ion . 

The  above  protocol  between  ainole-particle  solution  and  Faxen  law  is  dictated  by  the 
literature  —  single-particle  solutions  are  more  readily  available,  the  Faxen  laws  are 
not.  However,  for  new  shapes,  the  efficient  procedure  is  to  construct  the  appropriate 
tensor  Identity,  which  then  leads  to  both  the  single-particle  aolution  and  Faxen  law. 

For  prolate  spheroids,  the  singularity  solutions  of  Chwang  and  Wu  (1974,  1975)  furnish 
the  identities  necessary  for  the  derivation  of  the  Faxen  laws.  In  order  to  simplify  the 
algebraic  manipulations,  their  various  cases  have  been  combined  by  using  dyadic  notation, 
into  three  solutions  for  a  stationary  spheroid  in  a  uniform  stream,  vorticity  field  and 
rate-of-strain  field.  Wie  conversion  to  dyadic  notation  is  straightforward  for  the  uniform 
stream  and  vorticity  solutions,  since  the  solutions  correspond  to  cases  where  the  spheroid 
axis  jJ  is  either  parallel  or  perpendicular  to  the  relevant  vector  (the  uniform  stream  or 
the  pseudo-vector  for  the  vorticity  field).  For  the  rate  of  strain  field,  the  three  cases 
correspond  to  axisynnetric  straining  and  hyperbolic  straining  in  planes  containing  and  per- 
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pendicular  to  the  spheroid  axis*  The  tensorial  form  can  be  obtained  either  from  geometri¬ 
cal  arguments,  by  inspection  (Batchelor  and  Green  1972a)  or  after  rigorous  algebraic  steps 
(Brenner  1974). 

The  solution  of  Chwang  and  Wu  for  a  stationary  prolate  spheroid  (major  semiaxis  a, 
focal  length  2c,  eccentricity  e  =  c/a,  and  axis  d)  are  shown  below.  The  constants 
are  given  in  Table  1. 


Table  1.  Constants  for  the  velocity  representation  for  the  spheroid 
i)  Constants  derived  from  Chwang  and  Wu  (1974,  1975) 

a1  *  e2|-2e  +  (1+e2)  log(^~ 

a2  -  2e2|2e  +  <3e2-1)  logfj^)! 

y  =  ( 1-e2  ){  2e-(  1-e2  )  loa(|i|)| 

Yj  »  (1-e2){-2e+(1+e2)  iog(-i^S)) 
y’3  =  Y3(1-e2)-1 

a3  -  2e2Y3l-2e  +  logf  ~)  }  j  2e  ( 2e2 -3 )  +  3(1-e2)  log(~)} 
aj  =  e2Y3 1 -2e+(  1-e2  )  log(-|-^')  j  (  2e  ( 2e2 -3 )  +  3  (1-e2)  l°g(-3~)} 
a4  =  2e2  ( 1-e2  )  f  2e(3-5e2  )  -  3<1-e2)2  log(-J— )} 
a5  =  e2  j 6e-  (  3-e2)  log(^~]}-1 


ii)  Constants  used  in  equations  (3.7a-c),  (3.8a-c) 
a'  =  a3  +  »  e2( -2e+ ( 1 +e2  )  l°g(y^)) 

Y'  5  Y3  +  Y3  »  (2-e2)!-2e+(1+e2)  l°g(-|— 
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v(x)  «  U*  -  U* • (a  .dd  +  a„(6  -  dd))  . 

^  ^  ^  1'^'  2  *  ^ 

SI  (l  +  (C2  -  i2)  v2)  I(x-£)d5 

4e 


(3.7b) 


+  exud'd  x  /f  ( +  (c2-f2)  V2}d*e~(£)d5 

C  8e  ~ 

-  ~  xua3e3{ydd  +  Y'(S  -  dd)}*M 

Sii  -  8*"f-  5  Wi  -  jW* '  i4k*> 

- 1  ♦  di*jtdk  ♦  «it«jdk  +  sik«y*  ~  4didjdkdt> 

'  I  Wjt  +  4i«4jk  -  4ii4kl  +  fidi4kl  +  4ijdkdl  +  VjVt 

-  di4  jtW jkV4ikd  jV4 itd jdk 1  }^-c  (=2-52){1+(c2-€2) 


(3.8b) 


2~*2'  ill®!  V2)e“  (£)d? 
8e 


(3.8c) 


-  W(d.e  d,  ♦  ^ikl^i:c<c2<2){^  x  v“<£)}kd?  . 


These  expressions  for  the  Faxen  law  involve  integrals  of  the  ambient  velocity  field 
and  its  lower  order  derivatives,  whereas  the  ones  developed  by  Brenner  (1964b)  and  Rallison 
(1978)  use  infinite  series  involving  derivatives  of  the  ambient  velocity  field.  If  the 
ambient  velocity  field  is  derived  from  a  numerical  method,  the  latter  approach  may  run  into 
problems  caused  by  the  loss  of  numerical  resolution  (finite  difference  solutions)  or 
differentiability  (finite  element  solutions).  Our  integral  forms  yield  the  known 
differential  forms  if  one  expands  v,  VX.y  and  g  in  Taylor  series  about  the  spheroid 
center. 

For  the  important  problem  of  determing  the  motion  of  a  force-free  and  torgue-free 
spheroid,  (3.8a-c)  shows  that  the  translational  and  rotational  velocities  are: 

(3.9a)  0  =  1-  J®c{l  +  (c2-?2)  ( -  V2}v”(£)dC 

4e 


3  ,c  2  2  •» 

ui  =  ~  C  (c  -5  >  V  x  v  (UdS 

~  3  -C  ~  ~ 


(3.9b) 


8c 


+  -L.  ■■  .e..„  (c2-C2)(l  +  (c2-C2)  V2)d  x  e"(£,)*d  d f. 


3  2  ;  -c 

4c  (2-e  ) 


8e 


From  (3.9b),  the  change  in  the  orientation  of  the  spheroid  follows  as: 


-1 1- 


(3.9c) 


3  lc  <c2-f2)nV)  «  a  ac 

*  -e  ~  ~  ~ 


4c 


3  ;  -c 


4c  { 2-e  ) 


■z~  Sc  (c2-52){1  +  (C2-?2) 
2  v  ~c 


-■  72}  (g“(C  )«d  ~  e"(g  );ddd)dC  . 

8e  ‘  '  ~ 


For  homogeneous  ambient  velocity  fields  (constant  velocity  gradient)  the  vorticity  and 
rate-of-strain  tensors  can  be  taken  out  of  the  integrals  and  (3.9b)  and  (3.9c)  reduce  to 
the  Jeffery  equations  as  given  by  Leal  and  Hinch  (1972).  The  application  of  equations 
(3.8a-c)  and  (3.9a,b)  in  the  method  of  reflections  (the  module  integration  step)  will  be 
the  subject  of  parts  II  a^id  IV. 

There  are  two  complications  in  the  generalization  of  the  procedure  to  particles  of 
arbitrary  shape.  Because  of  the  high  degree  of  symmetry  possessed  by  spheres  and 
spheroids,  relatively  few  solutions  (or  tensor  identities)  were  required.  The  spheroid 
involved  two  solutions  for  the  uniform  stream,  two  for  the  vorticity  field  and  three  for 
the  rate-of-strain  respectively.  For  the  general  case,  three  streaming  solutions,  three 
vorticity  solutions  and  five  rate-of-strain  solutions  are  required.  Secondly,  it  is 
unlikely  that  analytical  expressions  will  be  available  for  these  eleven  single-particle 
solutions  (or  identities). 

The  solution  technique  of  Youngren  and  Acrivos  (1975)  is  recommended  since  it  will 
lead  to  Faxen  laws  involving  integrals  of  the  ambient  velocity  field  over  the  surface  of 
the  particle.  They  solved  the  integral  equation  for  the  velocity  field  by  collocating  the 
surface  of  the  particle  into  small  subdivisions.  The  approximants  for  the  velocity  in  the 
uniform  stream,  vorticity  field  and  rate  of  strain  field  cases, 

v ( x )  »  (T  -  ~  L  f(x' CKx-x*  )dS(x’ ) 

~  ~  ~  8*li  rS  ~  ~  m  ~  ~  ~ 

(3.10a) 

with  f(j')  *  U  *£(*')  • 

v<u>  -  S"  ■  *  -  is?  i<5’>-i«5-*,)<3s(x', 

(3.10b) 

with  fix'  )  -  0“ *C(x' )  , 
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(3.10c) 


v(x)  =  E*x  -  — -  f  (  x1  )•  I  (x-x*  )dS(  x'  ) 

~  ~  —  ~  8tt  y  S  ~  ~  ~ 

with  ,£(£'  )  =  )  , 


where  the  resistance  tensors  C  and  M  depend  on  the  particle  shape,  lead  to 

"approximate"  Faxen  laws  for  the  force,  torque  and  stresslet. 

(3.11a)  F  =*  v"(x'  )*A(x'  )dS  -  V‘j>s  A(x'  >dS  -  (w  *  x'  )*A(x'  )dS 

(3.11b)  T  *  v  (x*  )*C(  x'  )dS  -  U*£g  CUMdS  ~  f>s  (u  x  x'  )*C(x'  )dS 

(3.11c)  |  -  v*"(x*  )«S<x'  )dS  -  U»£g  H(x*  )dS  -  £g  (u  x  x'  )*B(x‘  )dS  . 


4.  MODULF  INTEGRATION:  HI  BETWEEN  TWO  STATIONARY  SPHERES 

The  method  of  reflections  and  the  Faxen  laws  of  the  previous  sections  are  combined  in 
this  section  to  calculate  the  drag  on  two  spheres  of  radii  a  and  b,  with  centers  at 
and  jjz,  held  fixed  in  a  uniform  stream  U  .  The  drag  on  sphere  1  as  calculated  by- 
standard  reflection  techniques  in  Happel  and  Brenner  (1965)  is: 


(4.1) 


El 

6  irya 


3  t  *  <L  2b  _  (!  x  ,  27  a2  +  1  3)  ^ 

4  R  16  r2  V4  64  4  R3 


81 

256 


-  dd> 
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3  b 
2  R 


9  ab 


4  R2 


&  X2 
8 


il  x2 

16 


~!dd|*u" 


R 


where  d  «  (lS2”X , )/ l22"-5l  *  •  *  “  b/a  an<1  R  =  !2J2--Si  *  *  The  method  presented  here  not  only 
recovers  this  result,  but  solves  the  parallel  and  perpendicular  cases  simultaneously. 

The  zero-th  reflection  at  each  sphere  yields  the  Stokes  drag  as  the  leading  contribu¬ 
tion  on  each  sphere: 


e;0’ 


(4.2a) 


6wyaU 


T1 


(4.2b) 


(0 )  • 
Z2  “  6*WbU 


The  reflected  fields  are  the  Stokes  disturbance  fields > 

(4.3a) 


-  a2  .2  £(2“2l ’ 

v.  -  -6»uaU  *(1  ♦  ~~  V  ) 

~1  ~  6  8*|i 


(4.3b) 


.  b2  2  1<S-Z 2) 

v  -  -6wybU  •  (!  +  §-  V  > 


Buy 


.<»> 


The  contributions  from  the  first  reflection  at  sphere  1,  F  ,  is  obtained  by  applying 
the  Faxen  law  (for  the  force)  with  v2  as  the  incident  field. 


F?1’  -  6xya(  1  +  V2)|-6»ub()"«(1  +  r- 

"'l  O  o 


S2  „2 


b2  „2 .  J(*"*2  ,I 


Sty  lx-*. 


so  that 
(4.4) 


2  .2. 


rn)  .  . 

ll _ 2v„*.(i  *  <•  3  *2 

6*Ua 


*-jbU*{1  + 


(V  I  •  0  since  J  is  a  solution  of  the  Stokes  equations).  The  substitutions  for  I  and 
V2X  reduce  equation  (4.4)  to 


6tya 


(4.5) 


> 

♦  [-f|+  (jX  ♦  ix3)]dd|.u“  . 


(2) 

The  contribution  fro*  the  second  reflection  at  sphere  1,  F  ,  is  obtained  by 
applying  the  Faxen  law  with  y,2  (which  was  created  by  reflecting  1  at  sphere  2)  as 
the  Incident  field.  Thus  the  multipole  expansion  for  v12  ia  required  accurate  to 
0(F'4): 

v  .  jn  ll  +  ,s< 1 )  +  Tn>,  lUffcjjgal  .  Qn>  ■U>,.M,(~"~2.> 

12i  2j  8»|i  2jk  2jk  8»u  2jkt  16*u 

and  the  contribution  from  the  second  reflection  follows  as 
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.(2) 


2i 


(4.6) 


f  1 )  »2  2  ^ 

- -  -  ’  (1  +  | ~  V2)  - =-| 

6nya  23  6  8iru  |x=x, 


+  (sn>  +  Tn> 

2  jk  2  jk  8k|i  |x=x, 

-  o'1*  111,kt<iS~~2) 


2jkt 


16k  y 


+  0(R_6) 


(1) 


The  momenta  in  (4.6)  must  be  simplified  to  obtain  the  final  expression.  F^  is 
obtained  by  switching  a  and  b  in  (4.5).  After  the  substitution  for  J  and  F*1**  t*'e 
first  term  in  (4.6)  reduces  to 

(4.7a)  {(fg  ^1+  (|  X  t  ^  X3)  ^](«-dd)  +  [|2|  -  (|  V  +  |  X3)  ^ldd[.o“  . 

R  R  R  R 

In  the  second  term  of  (4.6),  the  dipole  moments,  §  and  £  are  needed  accurate  only  to 

The  Faxen 


—2  —2 
0(R  )  since  they  are  multiplied  by  the  Stokes  dipole  which  decays  as  R 


laws  for  these  moments  yield: 


.(1)  1l-i>k(*~Z2> 


2  jk 


8i»u 


20  _._v3  1lj,k(~2~'g1> 


3  ’Mbeijk 


8s  u 


(1)  Iij.k(^~»2) 


3  3ii  k^2*"^1^ 
T'  '  - —I  -  -4syb  B  -ilti — - 

*2  jk  8su  1,5=8,  Ijk  8sU 


where  e,jk  =  j  (v  ^  +  v.^  ^)  and  0  =  ^  ^v1k,j  “  v1j,k^  are  the  rate  of  strain  and 

vorticity  fields  of  the  Stokes  solution.  The  substitutions  for  e,,  0,  and  the  Stokes 
dipole,  simplify  these  contributions  as 


(4.7b) 


X3  dd*u" 

4  4  ~ 

R 


3-  x3  («-dd)*u“ 

4  r4  - - 


and 
(4.7c) 
respectively. 

Hie  last  term  in  (4.6)  makes  a  contribution  through  an  0(R-^)  term  in  the  quadrupole 
moment  Q: 
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The  Faxen  law  for  Q  ia  shown  In  Kim  <1983)  along  with  a  recursive  formula  for  the  Faxen 


law  for  the  n-th  multipole  moment-  However,  the  leading  order  term  in  Q  can  be  deduced 

( 1 )  2 

from  its  presence  in  the  Stokes  solution.  After  the  substitutions  tor  F 2  and  7  J  the 


last  contribution  simplifies  to 
(4.7d) 


<16  *3  (4tK>  -  I  a3 

R  R 


The  calculation  of  the  contributions  from  the  third  and  fourth  reflection  is  easier 
since  only  the  leading  order  terms  are  needed  at  each  reflection  (point  force 

m 

approximation).  As  the  reflections  are  traced  back  to  the  ambient  field  V  ,  each 
reflection  contributes  a  factor  of  0(R-1 )  from  the  decay  in  I ( x ^  — x 2 ) . 

-<3)  2 

(-  7  b)  (-  i 


6wya 


(4.8) 


»"  S  >  <*-dd> 


27  ,2  .  - 

X  ■  dd)*p 
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,<♦> 


&-(-ib)(-J.)  I4tx2-*1)-0" 


(4.9) 


(81 


256  *2  $  ♦  S  *2  $  af 

R  R 


The  contributions  from  the  sero-th  through  the  fourth  reflections  given  in  (4.2), 
(4.5),  (4.7a-d),  (4.10)  and  (4.11)  sum  to  (4.1).  As  mentioned  earlier,  the  method 

m 

generates  simultaneously  the  solutions  for  streams  parallel  (dd*  U  terms)  and  perpendicu¬ 
lar  ((5  -  dd)*0  terms)  to  the  sphere-sphere  axis.  The  method  introduces  even  greater 
savings  in  the  mobility  problems  of  the  following  section.  The  parallel  and  perpendicular 
cases  will  be  handled  simultaneously,  but  more  importantly,  the  method  solves  directly  a 
specific  mobility  problem.  For  example,  the  motion  of  torque-free  particles  (e.g., 
sedimentation)  will  be  solved  by  imposing  the  torque-free  condition  at  each  reflection. 


The  result  will  be  the  one  of  interest,  and  thus  other  subsidiary  (and  unwanted)  problems 
will  be  avoided. 


5.  MODULE  INTEGRATION:  SEDIMENTATION  OF  TWO  SPHERES 

The  problem  of  determing  the  motion  of  two  torque-free  spheres  under  the  influence  of 
external  forces  (e.g.,  sediraenation)  can  be  considered  the  inverse  of  the  problem  discussed 
in  the  previous  section  (Brenner  1964a).  Our  immediate  concern  here  is  to  show  that  the 
present  method  gives  results  that  agree  with  the  expressions  for  the  mobility  functions  as 
deduced  from  the  expressions  given  by  Happel  and  Brenner  (1965). 

Consider  two  spheres  and  the  geometry  of  the  previous  section.  External  forces 
and  but  no  torques  are  imposed  on  each  sphere.  The  objective  is  to  express  the  trans¬ 

lational  and  rotational  velocities  U.,  U  and  iu<(  u,  in  terms  of  F,,  F„  and  d. 

~  1  *2  1  **  I  ~ 

Since  this  form  is  not  the  one  given  in  Happel  and  Brenner,  theirs  will  be  rearranged.  We 
start  with  the  general  relation  between  the  forces,  torques,  translational  velocity  and 
rotational  velocity  given  by  Brenner  (1964a): 


$11 

$12 

£11 

£12 

$21 

$22 

£21 

£22 

£n 

£12 

2n 

S12 

£21 

£22 

§21 

§22 

TOe  matrix  elements  are  second  order  tensors  and  the  usual  rules  for  matrix  operations  are 
in  effect.  For  the  problem  of  interest,  -  Tj  -  0,  so  that  and  may  be 

eliminated  from  (5.1)  to  give: 

r_  _.-i1  r. .  . . r..'i 


Ff  <6*ua) 


£2<6*ub) 


A*  A* 
*11  *12 


K.  A’ 
-21  -22 


2i  “  22  “  0  • 


The  symmetry  in  the  two-sphere  geometry  leads  to  a  further  simplification  --  each 
^g(a,8  -  1  or  2)  can  he  written  in  terms  of  scalar  functions  for  the  parallel  and 
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perpendicular  cases,  i.e 


(5.3)  ^  -  Xag(a/F,X)£d  ♦  *aB<a/*»x><*-jM)  . 

Since  the  tensors  coauiute,  the  inverse  of  the  matrix  in  (5.2)  can  be  written  as 

-1 


(5.4) 


&22-B 

•a?2 

-AI.'D*1 

A* 

L  "21  - 

■n 

-1 


with  0  *  fi}i"£22  ”  »2 1  * *1 2 '  We  now  substitute  Happel  and  Brenner's  (1965)  expressions 
(6-3.51)  and  (6-3.96)  for  X  and  Yi 
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(5.5a)  X 


11 


i  +  H*  (-f  »a  V  - 

R 


3  ,2  .  9  ,3,  r 
R4 


(5.5b)  X 


.  .  I  &  +  fl  x  _  27  2  I  3^  _  9  f  2  27  3  a 

12  2  R  l2  6  2  >  3  4'*  *8  *  *  K  > 


♦  .. . 


(5.5c)  Y„  -  1  ♦  fe  +  I  +  32  ~4  +  “* 

R  R 

„  3  b  1  /.  .  27  2  3^  a^  27  ,  2  9^  3  41  af 

(5.5d)  »12  ’  -  4  i  *  4  +  X  +XJj"64(X  +  16  X  +*J  • 

R  R 

The  expressions  for  X22,  X21,  Yj2  and  »21  are  obtained  by  switching  a  and  b  in 
(S.5a-d).  The  matrix  elements  in  (5.4)  can  now  be  written  as 

.-1, 


(5.6) 
and 

(5.7) 


+  X22<X11X22  -  X21X12>_1^  ♦  Y22(Y11Y22  ’  Y21Y12>" 


-ai2-S_1  ’  -X12<X11X22  *  X21X12,_1^  "  Y12<Y,,Y22  '  Y2 1Y 12 1 '  ’  » 


so  that  after  some  tedious  algebra,  the  translational  velocity  can  be  written  ass 
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H,  -  It1  +  ^2  ~  4^  x3  4  +  •••153  +  f1  +  ^2  +  ^4  +  —1  <«■«>}• 


0_ 

R"  R4 
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6*ua 
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♦  [j  i + j  <*+*3>  4  ♦  %  ♦  •••]<«-«>}•  12 

R  R 


6k  ub 


The  zeros  that  appear  in  (5.8)  are  not  accidental,  but  are  consequences  of  the  reflection 


process  and  the  properties  of  the  Oseen  tensor,  as  shown  in  the  following  direct 
computation  for  U ^ . 

We  now  obtain  (5.8)  with  the  present  method  by  writing  as  the  sum  of  the 

contributions  from  all  reflections  as  obtained  from  the  Faxen  law: 

S;°>  -  Foetus)-1 


si" 


(l  +  ~~  V2)v,|. 

v  6  J~2  I  X"X1 


si21 -(f  a- 


a2  „2s 


(1) 


Since  y2  is  the  Stokes  solution,  U 1  simplifies  as 


a2  „2 


.2  ,  I  ( x-x ) 

b  «  ~  ~2 


8»u  |  x-x 


(5.9) 
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6syb 


2 . .  2. 


(3  b)(l  + 


•?Sb  •  l  - 1 <x+x3>  -  f-j  •$ *  i  <>+*3> 


(2 ) 

The  calculation  for  U.j  introduces  a  novel  twist  which  turns  out  to  be  a  labor  saving 
feature  in  many  mobility  problems.  Since  the  first  reflection  already  contributed  p^, 
subsequent  reflections  cannot  contain  a  monopole  term.  Furthermore,  since  the  spheres  are 
torque-free,  the  multipole  expansion  for  any  of  the  higher  order  reflection  fields  must 
lead  off  with  the  stresslet  term: 


(5.10) 

with 


v.2  •:  xu 

S  (i  +  ^  v2)  - AJxL +  ... 

jkv  10  ’ 


8»U 


When  such  reflection  fields  are  used  as  incident  fields  at  the  next  reflection,  the  leading 
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mi. 


order  tan*  in  the  translational  velocity  at  that  next  reflection  will  be  of  0(R"3) 


•nailer  than  the  leading  term  in  the  tranalational  velocity  contributed  by  the  previoue 
reflection,  i.e.,  U<1n+1  -  0(R-3 Ig*"1 .  The  decay  in  the  Stokea  dipole  contribute*  a  fac 

of  o(R-2)  and  the  *tre**let  i*  OfR-1)  smaller  than  the  tranelational  velocity  of  the 

are  of 


preceding  reflection.  Therefore,  the  leading  order  term*  in  U®2'  and  uj  * 


0(R-4)  and  0(R-7)  respectively. 


(5.11) 
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Ix-x  1 
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The  translational  velocities  from  these  contributions  sum  to  (5.8)  and  as  mentioned 
earlier,  the  *ero  coefficients  in  the  mobility  functions  are  due  to  the  weaker  interactions 
beyond  the  first  reflection. 

The  rotational  velocity  jjj1  can  be  related  to  dX£1  by  inserting  (5.8)  into 
expressions  (6-3.97)  of  Happel  and  Brenner  (1965), 


(5.13) 
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The  result  1st 


(5.14) 
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again,  the  sero  coefficients  are  due  to  the  weaker  interactions. 

Our  method  of  reflections  sums  the  contributions  from  each  reflection.  Again,  it 
would  appear  that  factors  of  0(R-3)  will  be  introduced.  However,  the  stresslet  field  is 
irrotational  so  that  the  Interactions  are  even  weaker  (at  most  0(R-5)).  Therefore, 
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(5.15a) 


-£<y.v2> 


=  d*F  (6itgb)  1 

~  ~2  4  r2 


(5.15b) 


7i < x-x0 ) 

w*2*  .1  (Vxv  }J  =  -  ?x|s  1  : - 2—  * 

~1  2  ~12  I  x=x  n  2  =2  Situ 


d»F  (6wua)_1[-~  +  .. 
1  R5 
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~1  2  ~212  |x=x1  2  '=2  8nu  1 


=  dxF2(6*ub)_1  [^g 
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and  the  rotational  velocity  follows  as 


-It  0 

*F,(6llga)  —  + 
!S,  '-d~1  V 


(5.16) 
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In  summary,  the  method  reproduces  the  resistance  and  mobility  functions  for  spheres 
given  by  Happel  and  Brenner  (1965).  Furthermore,  the  results  are  consistent  with  the 
recent  and  more  extensive  calculations  of  Jeffrey  and  Onishi  (19S3).  In  the  following 
parts,  the  method  will  be  applied  to  HI  problems  involvina  spheroid-spheroid  inter¬ 
actions,  spheroid-wall  interactions  and  multi-particle  (three-sphere)  interactions.  These 
calculations  will  lay  the  foundation  for  future  investigations  into  the  effect  of  multi- 
particle  HI  on  the  rheology  of  a  suspension  of  nonspherical  particles. 
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NOTATION 


a  radius  of  sphere  at  x^. 

4  element  of  the  resistance  tensor, 

b  radius  of  sphere  at  X2> 

B  element  of  the  resistance  tensor, 

c  distance  from  center  to  foci. 

C  element  of  the  resistance  tensor. 

d  unit  vector  denoting  orientation  of  spheroid  axis. 
e  eccentricity  of  the  spheroid. 

£  rate  of  strain  tensor. 

£  force  exerted  on  the  particle  by  the  fluid, 

i  Oseen  tensor  defined  hy  equation  (2.2). 

U  outward  normal  vector  for  a  surface, 

p  pressure. 

quadrupole  moment  of  the  surface-force  distribution. 

R  center  to  center  separation  between  two  spheres. 

5  stresslet  or  symmetric  part  of  the  stress-dipole. 

T  torque  exerted  on  the  particle  by  the  fluid. 

0  particle  translational  velocity, 

v  velocity, 

x  position  vector. 

jf'  point  on  the  surface  S(x'). 

X  scalar  function  in  the  resistance  tensor,  parallel  problem. 

Y  scalar  function  in  the  resistance  tensor,  perpendicular  problem. 

a  constants  in  the  Chwang-Wu  singularity  solutions. 

Y  constants  in  the  Chwang-Wu  singularity  solutions. 

S  identity  tensor. 

E  alternating  tensor. 


-23- 


liiisaimiirnrii  i  rr  'ii'lmn  rutUm 


X  ratio  of  sphere  radii,  b/a. 
u  viacoaity. 

C  vactor  da noting  poaition  on  the  spheroid  axis. 

2  stress  tensor, 

w  particle  angular  velocity. 

2  vorticity. 

Q  vorticity  tensor. 

Subscripts 

1.2 

i. j.k.t.m 

Superscripts 

(n)  denotes  association  with  the  n-th  reflection. 


refers  to  spheres  at  ^2* 

indices  used  in  the  Einstein  summation  convention 


ambient  field 
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